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Abstract 

We generalize the concept of optical Galton board (0GB), first proposed by Bouwmeester et al. [Pliys. Rev. A 61, 013410 
(2000)], by introducing the possibility of nonlinear self-phase modulation on the wavefunction during the walker evolution. 
If the original Galton board illustrates classical diffusion, the 0GB, which can be understood as a grid of Landau-Zener 
crossings, illustrates the influence of interference on diffusion, and is closely connected with the quantum walk. Our nonlinear 
generalization of the 0GB shows new phenomena, the most striking of which is the formation of non- dispersive pulses in the 
field distribution (soliton-like structures). These exhibit a variety of dynamical behaviors, including ballistic motion, dynamical 
localization, non-elastic collisions and chaotic behavior, in the sense that the dynamics is very sensitive to the nonlinearity 
strength. 

PACS numbers: 03.67.Lx, 05.40.Fb, 05.45.Yv, 42.65.-k 



On 

> 

00 
O 

o 
\o 
o 

-(— > 



X 



INTRODUCTION. 

The Galton board, or quincunx, is a matrix of regu- 
larly spaced pegs fixed to a board through which pellets 
fall impulsed by gravity. The final distribution of pel- 
lets' locations at the bottom of the device follows the 
binomial distribution, and thus the Galton board consti- 
tutes a realization of the random walk. The importance 
of random walks does not need to be emphasized here, 
as their presence is ubiquitous in science. They are im- 
portant, in particular, as a tool in classical computation 
(the best known algorithms for solving some particular 
problems are based on their use [1]). For sure, this is 
one of the main reasons behind the present interest on 
the quantum counterpart of random walks, the so-called 
quantum random walks [2|| or, more appropriately, quan- 
tum walks (QWs). Moreover, from a fundamental point 
of view the study of quantum counterparts of important 
classical phenomena, and viceversa, is of obvious interest. 

The QW has been introduced from several different 
perspectives: In the seminal papers (in 1993 Aharonov et 
al. [2| introduced the QW as a generalization of the ran- 
dom walk, and in 1996 Meyer [3] introduced it as a non- 
trivial quantum cellular automaton) the computational 
aspects were not stressed, but later Watrous [4] inde- 
pendently introduced QWs from a quantum algorithmic 
point of view, and other versions of the QW (the so- 
called continuous QW, that can be viewed as a quantum 
generalization of the Markov chain) were also proposed 
[5|. Today there is a considerable amount of papers de- 
voted to QWs, and we refer the reader to existing reviews 

Not only one can think of quantum versions of the 
random walk, one can also think of wave [8|| as well as 
quantum versions of the Galton board [9^. The differ- 
ence between the wave and quantum versions lies in that, 
in the wave version, it is classical waves (e.g. optical 
waves) what are used, while the quantum version uses 



amplitude probability waves (wavefunctions) . The so- 
called optical Galton board (0GB) was first proposed by 
Bouwmeester et al. [8], and was introduced as a grid 
of (optical) Landau-Zener crossings. Bouwmeester et al. 
showed, both theoretically and experimentally, the exis- 
tence of spectral diffusion, as well as dynamical localiza- 
tion in their particular proposal for an 0GB. As for the 
quantum version, the quantum quincunx of Ref. [9] is a 
quantum-optical proposal for the implementation of the 
QW. 

Although classical waves and wavefunctions are differ- 
ent in a deep sense, a very interesting point raised by the 
0GB is that it can be understood, to some extent, as an 
optical realization of the QW [lO, UjJ, ll2|- It is conve- 
nient to stress that there are small differences between 
the 0GB of [8] and the QW, but as it was shown in [lo|, 
the 0GB reduces to a QW with an appropriate param- 
eter setting of the device. Moreover, Wojcik et al. [l2| 
suggested that their generalization of the Q W (consisting 
in the introduction of some additional position- depen- 
dent phase changes of the walker, see also [13|, [iJ, Il5| ) 
qualitatively describes the 0GB of (sl, as it reproduces 
the observed dynamical localization. These generaliza- 
tions of the QW have shown unsuspected connections of 
the QW with Anderson localization [l3| and quantum 
chaos [i2|,[l3. 

Here we propose a nonlinear generalization of the 0GB 
(NLOGB), consisting on the introduction of nonlinearity 
in the evolution of the walker. Given the connection be- 
tween the 0GB and the QW mentioned above, one could 
say that we are proposing a nonlinear generalization of 
the QW (the nonlinear QW). However, as we discuss 
below, our proposal makes full sense only from a clas- 
sical perspective, and thus our preference for the name 
NLOGB (nonlinear optical Galton board). As expected, 
the nonlinearity deeply modifies the QW dynamics, giv- 
ing rise to new and interesting phenomena which we in- 
vestigate in some detail. 
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After this introduction, the rest of the article is orga- 
nized as fohows: In Section II we briefly review the QW, 
as we win use its formahsm for introducing the NLOGB; 
in Section III we introduce the NLOGB as a nonhnear 
QW; in Section IV we describe the formation of sohtonic 
structures; in Section V we analyze the dynamics of the 
system describing the different phase transitions we have 
observed; and in Section VI we give our main conclusions. 



THE COINED QUANTUM WALK. 

Here we deal with the coined, discrete QW, in one di- 
mension. This process is better introduced as a quantum 
generalization of the random walk: In the random walk 
the walker moves to the right or to the left, depending 
on the output of a random process (e.g. tossing a coin); 
then the QW mimics the random walk in the existence 
of a conditional displacement that depends on the state 
of the coin, but differs from the QW in the fact that 
the coin is not a binary random process but a qubit. As 
the qubit can be in a superposition state, the walker can 
move simultaneously, say, in the two opposite directions. 
In order to make the dynamics nontrivial [3], the coin 
state must be changed (the analog of tossing the classi- 
cal coin) after each walk step, what is accomplished by 
the application of a suitable unitary transformation on 
the qubit. The main feature of the QW, as opposed to 
the random walk, is that the diffusion of theparticle is 
much faster (in the absence of decoherence [7]): While 
in the random walk the width of the probability distri- 
bution of the walker position grows as the square root of 
the number of steps, it grows linearly with the number 
of steps in the QW. Moreover, the probability distribu- 
tions have a very different shape (gaussian in the random 
walk, and resembling the Airy function in the QW). Let 
us now introduce formally the QW. 

The standard coined QW corresponds to the discrete 
time evolution, on a one-dimensional lattice, of a quan- 
tum system (the walker), coupled to a bidimensional 
system (the coin), under repeated application of a pair 
of discrete linear operators. Let Hw be the Hilbert 
space of the walker, with {|m) ,m G Z} a basis of Hw'j 
and let He be the Hilbert space of the coin, with basis 
{\u) ,\d)}. The state of the total system belongs to the 
space H = He ^ Hw cind, at time t, can be expressed as 
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The connection between states in consecutive times 
is made by an unitary linear evolution operator U, 
which can be written as U = UdUc, i-e., |^(0) = 
UdUc \ip{t — 1)). Here, Uc = C® I is the "coin toss oper- 
ator" with C G SU{2), typically chosen as the Hadamard 
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is the "conditional displacement operator", which moves 
the walker one position to the right or to the left, de- 
pending on whether the coin state is u or d, respectively. 
The main quantity related to the walk is the probabil- 
ity distribution function of the walker along the lattice, 
calculated as Pm{t) = \um,tf + \dm,tf = P^(t) + P^(t). 

We have already commented that the QW can be clas- 
sically simulated. In order to make things concrete, we 
consider the scheme depicted in Fig. 1, which repre- 
sents an optical cavity. A quasi-monochromatic light fleld 
enters the cavity through a partially reflecting mirror. 
When this fleld reaches the beam-splitter (BS in Fig. 1), 
it can follow two different paths, upper and lower in the 
flgure. These two paths play the role of the qubit (which, 
in this case, would be better called a cebit, following the 
terminology introduced in [17]), and the beam-splitter 
implements the unitary transformation [17] (the "coin 
toss" operator). Then in the lower (upper) path, the fleld 
frequency, which plays the role of the walker in this op- 
tical implementation, is increased (decreased) in a flxed 
amount Ao; by means of appropriately tuned electrooptic 
modulators. This is the flrst step of the QW. Then, the 
cavity mirrors reflect the light back to the beam-splitter 
and a new step of the QW is implemented, and so on and 
so forth. 

In this case, the QW occurs in the frequency distribu- 
tion of the output fleld, with the intensity of each fre- 
quency component playing the role of the probability of 
flnding the walker at a given position, i.e., P^{t) and 
P^(t) are spectral intensities in this classical-wave con- 
text, and not probabilities. In other words, after m cavity 
roundtrips, the spectrum of the output fleld exhibits the 
probability distribution of the QW. This is one of the 
schemes proposed in [l0| for the optical (classical) imple- 
mentation of the QW, where also the connection between 
the 0GB of Ref. [8] and the QW is given, and we refer 
the reader to that paper for full details on this type of 
classical (optical) implementation of the QW. Let us em- 
phasize that this scheme constitutes a realization of the 
optical Galton board. 

What this classical implementation of the QW (and 
others 18|, Ujl, l20|, l21|) suggests is that interference, and 
not entanglement, is the responsible of the QW charac- 
teristics. Entanglement would manifest in QWs in more 
than two dimensions, in the amount of classical resources 
needed for its implementation, as compared with a true 
quantum implementation, as already discussed in [10|. 
This does not mean that there is nothing quantum in 
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Figure 1: Optical cavity for the implementation of the 0GB 
and the NLOGB. EOMi and EOM2 are two electrooptic mod- 
ulators which are tuned for incrementing (decreasing) the field 
frequency in Acj. BS is a beam-splitter, and the cavity is con- 
stituted by four mirrors, one of which is partially reflecting 
and serves as input /output port. For implementing the 0GB, 
the upper and lower paths must be a linear optical medium, 
which must be replaced by a nonlinear optical medium (such 
as, e.g., an optical fiber) for implementing the NLOGB. 



the QW: It is the different physical meaning of Pm{t) 
(in a true quantum system, the probability distribution 
can be reconstructed only after a large enough number 
of measurements, while in the classical simulation the 
analog of the probability distribution corresponds to the 
field spectrum and can be seen completely at each walk 
step). The effect of decoherence could be different in 
classical and quantum implementations [7]. But, at least 
in the QW on the line, the quantum nature seems not to 
manifest, as it can be successfully simulated by classical 
means. See [22] for a discussion on this topic. 



INTRODUCING THE NONLINEAR OPTICAL 
GALTON BOARD. 



The optical cavity scheme of Fig. 1 serves us to in- 
troduce the nonlinear optical Galton board (NLOGB). 
It suffices to assume that light acquires some intensity- 
dependent phase while traveling through the upper and 
lower paths, i.e., that these paths are not made by a 
linear medium (vacuum), but with a nonlinear medium 
(e.g. a Kerr medium, like an optical fiber or similar). 
This is very easily taken into account with the QW for- 
malism introduced in the previous section that we will 
continue to use here: We only need to introduce one 
more operator describing the acquisition of the intensity- 
dependent (nonlinear) phase due to propagation in the 
Kerr medium, i.e., we have to generalize the unitary op- 



erator defined above in the following way: 

U{t) = UdUcUnl{t-l), (4) 

Uni{t) = Y. ^e'^^(™.*) \c,m) {c,m\ , (5) 

c=u,d m 

where Fc (m, t) (c = u^ d) is an arbitrary function of the 
probabilities (or intensities, in a classical context) P^{t) 
and P^{t) [23]. Notice that the role of Uni(t) is to add a 
nonlinear (probability dependent) phase to each of the 
spinor components. With the above formulation, the 
standard QW is obviously recovered when F^^ = F^^ = 0, 
and the generalized Q Ws of [l3| and 12|, llj| are recovered 
when Fu = Fd = m?(j)Q and Fu = Fd = m^o, respectively, 
with (j)o a constant phase. We see that a physical system 
like the one represented in Fig. 1 allows to implement 
a number of interesting generalizations of the QW in a 
relatively simple way. Let us emphasize that the 0GB 
of Bouwmeester et al. ^] is very close to what we are 
commenting [10|. 

In this article we shall consider one of the simplest 
forms for ([5]) by choosing Fc{m^t) = 27ra\cm,t\ (c = 
u^d), i.e., we assume that the nonlinear phase gained 
between two QW steps is due to a Kerr-type nonlinearity 
that acts separately on the two coin states {u and d) 
and has a strength a. The recursive evolution equations 
for the probability amplitudes can be easily derived from 
|V^(t + I)) = f7(t + 1) I V^(t)), yielding 
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As we show below, the nonlinearity just introduced 
deeply modifies the behavior of the probability distribu- 
tion Pm{t)- For this purpose, we perform a numerical 
study of Eqs. (|6|7p for different values of a. We shall 
consider a > for definiteness, since from Eqs. (|6)7p one 
easily sees that choosing a positive a, say a = aoj with 
some initial conditions {um,o]dm,o) is equivalent to tak- 
ing a = —ao and complex-conjugated initial conditions 
{um,o''> drn,o) • Morcovcr, we shall adopt, unless otherwise 
specified, symmetrical initial conditions localized at the 
origin, i.e., Um,o = ^mo/V^ and dm,o = iSmo/V^- 

From a classical (wave) viewpoint, the above process is 
a nonlinear optical Galton Board (NLOGB) and can be 
implemented with the same device we have commented in 
the previous section, provided that the two optical beams 
propagate in a Kerr medium (e.g., an optical fiber), as 
this nonlinear propagation exactly corresponds to what 
Uni represents. From a quantum viewpoint the imple- 
mentation of Uni is probably impossible because of the 



linearity of the Schrodinger equation. It is clear from now 
that the process we are proposing makes full sense only 
as a nonlinear 0GB, and will find conceptual difficulties 
as a nonlinear QW. 

In spite of the difficulties when speaking of a nonlin- 
ear QW, one should keep in mind that nonlinearities can 
be introduced in quantum systems through a clever use 
of measurement [2J, |25|, what keeps open the possibil- 
ity of implementing the proposed NLQW. Another, more 
realistic possibility concerns using systems described by 
nonlinear effective Hamiltonians, as Bose-Einstein con- 
densation, where QWs could be implemented [26|, or su- 
perconducting devices, just to mention a couple of poten- 
tial candidates. But these appear as remote possibilities, 
as compared with the immediacy of an optical implemen- 
tation in an optical device similar to that already used 
by Bouwmeester et al. [8]. 



FORMATION OF SOLITON-LIKE STRUCTURES. 

In Fig. 2 we represent the evolved probability distri- 
butions Pm{t) for a = (i.e., the standard QW) and 
a = 0.4. When a = 0, we observe the typical QW be- 
havior [6]: Pm{t) exhibits two peaks at the borders of 
the distribution, whose tails decay in the central zone, 
and whose maximum value monotonically decreases with 
time as the probability distribution broadens; and, most 
importantly, the width of Pm{t) is proportional to t. This 
probability distribution can be expressed, in some limit 
|ll| . as a combination of Airy functions propagating in 
opposite directions. 

The shape of Pm{t) for a = 0.4 is very different: Now 
the two peaks of Pm{t) contain most of the total proba- 
bility, around 30% each one in the case of Fig. 2, mostly 
distributed within a few lattice positions (see the inset 
in Fig. 2) . But the most striking characteristic of the 
probability peaks in this nonlinear case, is that their size 
and shape remain basically constant with time, except for 
small oscillations around a mean value. 

We will characterize the probability peaks by their 
position and intensity (i.e., the total probability they 
contain). As for the position, given the small ffuc- 
tuations on the shape of the peak, we use the "cen- 
ter of mass", defined as mcM = Zlrn^^^(^)' ^^^^ 
m e [mmax + Am,mmax - Am], mmax the position of 
the probability maximum and Am the width of the peak 
[27|. Only small quantitative differences are found be- 
tween the behavior of rrimax and that of mcM- 

The most important feature of the probability peaks 
is that they are non-spreading pulses, i.e., they prop- 
agate without distortion [28]. As these probability 
wave-packets do not spread in time, and present other 
particle-like features (see below) we can consider them 
as solitonic-like structures, and will simply refer to them 
as solitons. 
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Figure 2: (Color online) Probability distribution curves of 
Pm(t) for t — 300, with the initial condition iXm,o = ^mo/V2 
and dm,Q — i5mo/V^^ for a = (standard QW) and a — 0.4. 
The inset is a magnification of the right-moving probability 
soliton. Notice that Pm (t) is null for odd m (as t is even 
in this plot). We have represented only nonzero values and 
joined them for guiding the eye. 



Apparently, solitons do not require a minimum value of 
a to form: We have checked their existence for a > 0.01, 
and the analysis of the data from different (non-zero) val- 
ues of a does not suggest the existence of any threshold 
for the solitons formation. Nevertheless, the time needed 
for their formation (i.e., the transient until the inten- 
sity and shape of the probability structure is constant on 
the average) is larger for smaller a. This is appreciated 
in Fig. 3 (a), where the intensity of a single soliton is 
represented as a function of time for different values of 
the nonlinearity parameter a. Another important feature 
is that the width of the overall probability distribution 
Pm{t) or, equivalently, the soliton velocity, decreases as 
a increases, as shown in Fig. 3 (b), where the position 
of the solitons is represented for different values of a (in 
the case a = 0, where solitons do not exist, we have rep- 
resented the position of the center of mass of the max- 
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Figure 4: (Color online) Evolution of the ratio a/t for different 
values of a. 



Figure 3: (Color online) (a) Intensity (total probability) of 
the right-moving soliton. (b) Temporal evolution of the cen- 
ter of mass, mcM, of the right-moving soliton (the plot is 
symmetric for the left-moving soliton). The values of a are 
indicated in the plots. Initial conditions are as in Fig. 2. 



imum of Pm (t) for the sake of comparison). Therefore, 
solitons form after some transient, and are slower and 
more intense for larger a. This is the scenario we found 
for a < 0.474. 

In view of the phenomena described above, one might 
wonder whether the intrinsic quantum features of QWs 
are deteriorated or not, and if so, to what extent. One 
way to quantify the possible loss of the quantum ben- 
efits is by analyzing the time evolution of the standard 
quadratic deviation a = V< ra'^ > — < m > '^ . As al- 
ready discussed, the standard QW exhibits a character- 
istic a (X t. Given the transient which appears during 
the formation of the solitons, the question we ask our- 
selves is: does the quotient a/t go to a constant after the 
transient (i.e., for sufficiently large time), or will it decay 
slower? 

As can be seen from Fig. 4, the first possibility is in 
fact realized: after the transient, the standard deviation 
approaches the typical QW time evolution. Therefore, 



the long-term QW behavior is not degraded by the for- 
mation and propagation of the solitons. 



DYNAMICAL PHASES. 

We have been able to identify three different dynamical 
domains, or dynamical phases, in the behavior of solitons 
as a function of the value of a: Phase I, for a < aj c:^ 
0.474; phase II, for aj < a < an c:^ 0.6565; and phase 
III, for a > an. Let us describe these phases separately. 

In Phase I, the dynamics is very simple: Once solitons 
have formed, they exhibit the ballistic propagation al- 
ready shown in Fig. 3(b). Differently, in Phase II the two 
solitons start moving in opposite directions, as in Phase 
I, but after some time their velocity decrease till the soli- 
tons reach a turning point and then move backwards and 
collide at some later instant tcoi at m = 0. After the col- 
lision, the solitons continue moving apart indefinitely, as 
in phase I. An example of such behavior, for a = 0.49, is 
shown in Figure 5. Notice the appearance of small "com- 
munication packets" that are interchanged between the 
two solitons. Interestingly, the solitons intensity sharply 
decreases after the collision (for example, for a = 0.49 




Figure 5: (Color online) Color density plot showing the evolu- 
tion of Pm{t) as a function of t (horizontal axis) for a = 0.49. 
The vertical axis corresponds to the position on the lattice. 
Brighter regions indicate a higher probability. The two soli- 
tons are clearly visualized as intense strips. 
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the intensity of one soliton falls from 0.3062 before the 
collision, to 0.2426 afterwards), i.e., the collision of the 
two solitons is an inelastic one. Another feature that 
can be observed from the simulations is that, as a is in- 
creased from below aj (inside phase I), the intensity of 
the solitons increase up to a maximum value. It seems 
that the communication packets interchanged by the two 
solitons play the role of an attractive interaction, which 
is larger for larger intensities. This would explain the ex- 
istence of the above-mentioned turning point appearing 
at some critical value a/. Inside phase II, the solitons 
experience an inelastic scattering and loose a fraction of 
their intensity, which would prevent from recollapse. 

The method we used to determine this critical value, 
however, makes use of the fact that the collision instant 
tcoi decreases with a. Indeed, the function tcoi{(^) can 
be well reproduced numerically by a simple hyperbola 
^/tcoi = ci/oi-\-b (where the values of a and b are obtained 
by a numerical fit, with a coefficient of determination 
r^ = 0.99516, giving a = -0.0297 ± 0.0003 and b = 
0.0627 ± 0.0006). This numerically-obtained law allows 
to fix the frontier between phases I and II by the a value 
for which tcoi diverges (we obtained aj = 0.474 ± 0.007). 

As we made for phase I, it is worth investigating how 
the standard deviation evolves at long times, in order 
to quantify a possible departure from the characteristic 
quantum spreading. As before, we plot in Fig. 6 the 
quotient cr/t as a function of time for values of a corre- 
sponding to the second phase. Now the transient shows 
more complicated features, due to the recollapse of the 
two solitons (which manifests as the minimum appearing 
in both curves). However, as the solitons separate after 
the collision, the typical a (x t behavior shows up. 



Figure 6: (Color online) Same as Fig. 4, for values of a 
corresponding to phase II. 



Phase III, a > a//, differs from phase II in that, after 
the collision, the two emerging solitons do not necessarily 
collide or separate from each other. In fact, if a is in- 
creased beyond a//, the situation becomes quite compli- 
cated, as the evolution of the solitons becomes extremely 
sensitive to small variations in a. In this sense, we can 
say that phase III is a chaotic phase: For some values 
of a, the solitons become trapped and oscillate around 
the origin; with a slightly different value for a, however, 
the solitons eventually escape; and there are other a val- 
ues for which localization is found, which is characterized 
by an asymptotic setup of both solitonic structures at an 
equilibrium point. Interestingly, the latter possibility can 
occur at very distant site positions for slightly different 
values of a: For example, the right-moving soliton posi- 
tion oscillates between m = 5 and m = 9 for a = 0.6665; 
it remains static at position m = 162 for a = 0.6669; 
and again oscillates, around m = 7, for a = 0.6673. In 
Fig. 7 we show an example of the type of dynamics one 
encounters in phase III for the two values of a indicated 
in the figure caption. 

The results we have just described correspond to a par- 
ticular choice of the QW initial conditions, which guar- 
antees the symmetry of the probability distribution with 
respect to the starting position. In order to see how crit- 
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Figure 7: (Color online) Same as Fig. 5, but for a =0.6565 
(top) and a = 0.658197 (bottom). 



ical is the role of the initial condition, we have carried 
out numerical simulations for different sets of initial con- 
ditions, and have found that the dynamics is also very 
sensitive to this choice. Fig. 8 gives an idea of how dif- 
ferent things can be: we represent the evolution of the 
probability distribution for i^rn,o = ^m07<^m,o = and 
a = 0.2 (top) or a = 0.6 (bottom). For this initial condi- 
tion, the probability distribution is no longer symmetric 
(even in the standard Q W) , and this fact strongly affects 
the formation and dynamics of solitons. We are not go- 
ing to enter into an exhaustive description here; it will 
suffice to say that, in this case, there are also several 
dynamic phases: For small a, a single soliton forms, car- 
rying close to 60% of the probability, that moves like in 
Fig. 7 (top) (most of the rest of the probability is con- 
tained in small dispersive pulses that can be appreciated 
in the figure); for large a several solitons, with different 
intensities, can form, and localization phenomena similar 
to what we have described above can occur too, see Fig. 
8 (bottom). 



CONCLUSIONS. 

We have introduced a simple variation of the Opti- 
cal Galton Board (which can be understood as a clas- 
sical implementation of the discrete coined QW), based 
on the assumption that light propagates through a non- 
linear (Kerr-type) medium inside the optical cavity or, 
using the algebraic language of QW, based on the acqui- 
sition of non-linear -probability dependent- phases by the 
state during the walk. 




Figure 8: (Color online) Same as Fig. 5, but for a = 0.2 (top) 
and a = 0.6 bottom. The initial conditions are Um,o = ^mo 
and dm,o = 0. For a = 0.2 a single soliton is formed that 
carries 57% of the probability. For a = 0.6 this soliton has 
now a smaller intensity (32% of the probability) and becomes 
localized near m = 0, while a second soliton (20% of the 
probability) is formed. 



The most striking feature that the nonlinearity intro- 
duces, is the formation of soliton-like structures, which 
carry a constant fraction of the total intensity (probabil- 
ity) distribution within a non-dispersive pulse. We have 
characterized the dynamics of these solitons showing the 
existence of complex dynamics (from ballistic motion to 
dynamical localization) that is very sensitive to the ini- 
tial conditions. An important feature we found is that, 
in spite of the complicated behavior during the transient 
and possible recollapse of the solitons, the long term evo- 
lution still shows the characteristic QW feature in the 
cases when the solitons go away, in the sense that the 
standard deviation becomes a oct. 

It would be of greatest interest to have at hand an ana- 
lytical description of the solitons motion and interaction, 
specially during the formation transient and recollapse 
(when present), as done (approximately) in [10]. The 
additional complication due to non-linearities, however, 
makes this task cumbersome and lies beyond the scope 
of this paper. 

The described phenomena are, to the best of our 
knowledge, new in the field of quantum walks. The ex- 
citing features found here deserve, we believe, further 
research. 
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